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We have examined a Hubbard model on a chain of squares, which was proposed by Yajima et al as 
a model of an atomic quantum wire As/Si(100), to show that the flat-band ferromagnetism according 
to a kind of Mielke-Tasaki mechanism should be realized for an appropriate band filling in such a 
non-frustrated lattice. Reflecting the fact that the flat band is not a bottom one, the ferromagnetism 
vanishes, rather than intensified, as the Hubbard U is increased. The exact diagonalization method 
is used to show that the critical value of U is in a realistic range. We also discussed the robustness 
of the magnetism against the degradation of the flatness of the band. 

PACS numbers: 71.10.-w,71.10.Fd,75.10.Lp 



Recent progress in atom manipulation on solid sur- 
faces using scanning tunneling microscope has opened 
up an unprecedented possibilities for designed structures 
on atomic scales. Specifically, Hashizume et atJ has fab- 
ricated atomic wires from Ga adatoms bound to a row 
of dangling bonds on a hydrogen-terminated Si(100)2xl 
surface. Watanabe et afl have obtained its band struc- 
ture with the first-principles calculation. Interestingly, 
they have found that a partially flat band appears in the 
dispersion and suggested a possibility of the ferromag- 
netism arising from the flat band. Some of the present 
authorsa have then extended the electronic-structure cal- 
culation by going from a group III adatom (Ga) to a 
group V adatom (As). Unexpectedly, a flat (dispersion- 
less) band appears for some atomic configurations for the 
As wire as well, which again suggests a possibility of fer- 
romagnetism. As a possible origin of the flat band, Ya- 
jima et au have proposed a tight-binding model, which 
is thought to capture the essential feature of the arrange- 
ment of the arsenic atoms on the wire as a chain of con- 
nected squares, (fig.0) 

Possibility of ferromagnetism due to the electron- 
electron repulsion on flat one-electron bands, or a macro- 
scopic degeneracy in the one-electron states, was first 
shown rigorously for the Hubbard model by Lieba. While 
Lieb considered the models on bipartite lattices that have 
different numbers of sublattice sites CNa — Ye cx the 
magnetization), Mielketffl and TasakiQ have later con- 
sidered line-graph theoretical (non-bipartite) models in 
which flat bands appear from interferences between the 
nearest- neighbor transfer t and more distant transfer t' . 
The interference in each plaquette (which is triangular) 
is in fact a key to their rigorous reasoning why the flat 
band, when half-filled, leads to the ferromagnetism in 



the presence of the Hubbard repulsion U . One remark- 
able feature of the Mielke-Tasaki theorem is that an in- 
finitesimal U is enough to make the ground state fully 
spin-polarized despite the fact that the magnetism ap- 
pear from the electron correlation. 

Now, Yajima's model for the atomic wire is bipartite, 
so that Lieb's theorem is applicable if the site energies 
are all identical. Since the model is a chain of squares 
with Na = Nb, the ground state is spin-singlet at half- 
filling (number of electrons = number of sites). Hence, a 
kind of Mielke-Tasaki mechanism has to somehow work 
to realize ferromagnetism in this model. 

Usually Mielke's models (as exemplified by Kagomc 
lattice) and Tasaki's models are conceived as the graphs 
comprising triangles. Penc et au have actually shown 
that a chain of connected triangles is favorable for ferro- 
magnetism. One is then tempted to interpret the inter- 
ference mentioned above as an outcome of a frustration 
in odd-membered rings. So the challenge here is that: is 
it possible to have ferromagnetism a la Mielke-Tasaki for 
non-frustrated (e.g., a chain of even-membered rings as 
in the model proposed by Yajima et al). This is a dou- 
bly nontrivial question, since, even when one has a flat 
band from even-membered rings, it is by no means a suf- 
ficient condition for ferromagnetism and the applicability 
of the Mielke-Tasaki mechanism has to be examined. In 
this paper, we address to this question to conclude that, 
unexpectedly, a ferromagnetism a la Mielke-Tasaki does 
indeed exist for the non-frustrated model for the atomic- 
scale wire structure. 

We consider the Hamiltonian, which consists of the 
one-electron tight-binding part, Ho, and the Hubbard 
repulsion, Tiu, where 
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Here a ic ~ <ij CT create, respectively, electrons with spin 
cr at the four incquivalcnt positions a ~ d on the square 
(fig.|l|) in the ith unit cell, while t\ is transfer between Si, 
t 2 between As, and t 3 between Si and As. We have as- 
sumed a level offset, Ae between the As level and the Si 
level with n b i(T = b ia bi a etc being number operators. Here 
we adopt periodic boundary conditions for the system 
having L unit cells, and we take t\ = 1. 

This model has a special feature that, although the 
lattice contains only even-membered rings, we do have 
flat band(s) if the parametric relation, 

h = t 2 

t'3 = ti COS 9 

Ae = - 2sin6>) 

is fulfilled, where 9 is arbitrary (except for the condition 
9 / — 7r/2,0 as discussed later). It is in fact straightfor- 
ward to show that the single-electron Schrodinger equa- 
tion has L-fold degenerate solutions with the eigenener- 
gies 

£o = — 1 — sin 9. (1) 

Namely, this model has in general one flat band (the 
second from the bottom) for a given value of 9. We can 
alternatively place the flat band the third from the bot- 
tom (eq = 1 — sin9) with Ae = ti(l — 2 sin 9), which is 
equivalent to the former through gauge and electron-hole 
transformations. If one wishes to have Ae = 0, then one 
can set 9 = —tt/6, which is a special case of two fiat 
bands (at e = ±1/2) with an electron-hole symmetry. In 
figj^ we show typical band structures. The fact that the 
structure is constructed from even-membered rings has a 
reflection that, unlike Mielke's model or Tasaki's model, 
the flat band is not the bottom band, but a dispersive 
band lies below (unless one changes the sign of one of 
the transfers in the ring). 

Our first key finding is that the reasoning employed by 
Mielke and Tasaki in ref[9] can be applied to the present 
flat band. We can anticipate this, since the eigenfunc- 
tions on the flat band, 

= [C + sin ^4-cos0 4 

+ cos9a\ +la - sm9b\ +la - c| +1 J|0), 



have a notable property that each eigenfunction can- 
not be accommodated within a unit cell of the structure 
(a shaded area in fig. 1 ) . This feature is favorable for 
a realisation of the connectivity condition conceived by 
Mielketa 

We first apply the reasoning employed by Mielke and 
Tasakfa to the flat band here, while the effect of other 
bands are examined later. For this we can use the sub- 
space spanned by {(j>i,a} for the number of electrons fixed 
to L (i.e., the half-filled flat band). Any states in this 
subspace are represented as 

*= £ n # T n &i°> 

\L- t \ + \L l \=L i£LT jELl 

where L^, are arbitrary subsets of {1 ■ ■ ■ L} with 
and |LJ being the number of elements, respectively. 

If there is a state that satisfies TLu^gs = 0, then &cs 
is the ground state in this subscace and, following the 
argument by Mielke and TasakiH, we can show for any 
U > that such states exist, non-degenerate apart from 
the trivial spin degeneracy as follows. First we require 
that n%n°; l §GS = or nf^nf^os = for any i. To 
satisfy this condition, we must choose Lf and L| so that 
Lf PI ij, = 0. Hence, we can represent any ground state 
as 

*GS = E^)II4 W |0) 

a i 

where a represents the spin configuration. 

For Hu^gs — to be fulfilled we require n^n^$GS 
or n^nl^&Qs = for any i. This requirement amounts 
to 

where ^2 stands for a summation over a that satisfies 
a(i — 1) =| and <j(i) =|. From this we can see that 
g(v) = g(&x^ y ) for any x, y G {1 • • • L}, where a x ^ y de- 
notes the spin configuration obtained by exchanging d{x) 
and cr(y). The lowest energy state (Hu&gs = 0) with a 
given total S z (= \{L^ — L{)) is unique because all the 
g(cr)'s for each S z sector take the same value and fac- 
tor out. This is exactly the ferromagnetic ground states, 
since the fully polarized state does not experience U, and 
this concludes the proof of the above statement. 

In the present model a dispersive band lies below the 
flat one unlike the Mielke's or Tasaki's model. We can 
expect that a small enough U does not mix the bands 
since they are separated with a finite gap, but this has 
to be examined. For that purpose we take account of the 
one-electron states {ipi,a}i=i - L < on the dispersive band 
as well (where the number of states L < is just L here, 
since we have cply one dispersive band below the flat 
one). As Mielkeli3 discussed, sufficiently small U can be 
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treated with a degenerate perturbation theory, and we 
can show that the ground state of H = TLq + Jijj with 
L + 2L K electrons still has the total spin S — L/2 for a 
small U. 

By contrast, this approach becomes invalid for a large 
U, so that it is necessary to determine the total spin of the 
ground state from other methods. We have determined 
this numerically with the exact diagonalization. First, 
let us consider the case of Ae = 0. In fig.[| we show the 
result for the case of 6 electrons in 8 sites [L = 2) or 12 
electrons in 16 sites (L = 4), for which the band filling is 
n = 3/4 with the flat band being half- filled. We denote 
the lowest energy among the states having a given total 
S z by E min (S z ), and plot A = E min (S max ) - E m - m (0), 
where Sma* = L/2 is the S z when the flat band is fully 
spin polarized (Smax = 1 for the 8-site case, S'max = 2 
for the 16-site case), as a function of U. In agreement 
with the above analysis, the ground state has indeed the 
polarized flat band from an infinitesimal U . 

Our second new finding is that the ferromagnetism is 
destroyed when the interaction is too large (U > Uc)- 
The critical value Uc at which the ground state becomes 
unpolarized is seen to be Uc = 2.3 ~ 2.4, where the 
sample-size dependence is small. To the best of our 
knowledge this is the first example of the ferromagnetism 
that is unstable for large interactions. 

In order to check whether the transition to the un- 
polarized state is abrupt, we plot in the inset of fig.|| 
A = E min (2)-E min (l) along with A = £ min (2)-£ min (0) 
against U for the 16-site system. We can see that the 
transition from the state with S = L/2 directly into the 
5 = state. 

So far we have considered the case of half-filled flat 
band. In the As/Si system .this would require a doping 
with e.g., alkali metal atoms.Eil The ferromagnetism is ex- 
pected to be degraded when the flat band is pushed away 
from the half filling. It is believed that one-dimensional 
Tasaki's model is ferromagnetic only when the flat band 
is half-filled, while paramagnetic for lower electron densi- 
ties. In order to see if this is also the case in the present 
model, we calculate A = E min (l) — E m i n (0) as a function 
of U for the number of electrons decreased from 12 to 
10 (quarter-filled flat band) in the 16-site system. For 
the parameter region we have investigated, the ground 
state is paramagnetic, i.e., the model shares the insta- 
bility against the hole doping with the one dimensional 
Tasaki's model. 

We must also look into how the ferromagnetism pro- 
posed here is robust against the degradation of the flat- 
ness pof-the band. For Tasaki's model, Kusakabe and 
AokiL3'li3, and later Tasakio have shown that ferro- 
magnetic ground state survives small perturbations that 
make the flat band dispersive. Here we look at the sta- 
bility of the ferromagnetism against the perturbation, 

There are two reasons why we consider such an extra 



transfer. First, in real atomic wires, there should be 
a finite transfer between the adjacent As clusters, j-and 
secondly, the ferromagnetism on the ladder systemala is 
subject to recent investigations, so that the effect of a 
perturbation toward the ladder is of interest. In fig.^, we 
plot A = -Emin^max) — -Emin(0) as a function of U for 
St = 0.1,0.2 for 6 electrons in 8 sites. We can see that 
the ferromagnetism survives St = 0.1, while destroyed for 
St = 0.2. 

Next, we look at the robustness of the ferromagnetism 
for the general value of Ae / (9 / — tt/6). We can 
start from an observation that the ferromagnetism will 
be lost for 9 = or —tt/2, since {4>i.a} foils to satisfy 
the connectivity condition for 9 = invalidating Mielke- 
Tasaki's argument, or the lattice becomes disconnected 
with £3 = for 9 = —ir/2. We have numerically ob- 
tained Uc as a function of Ae for 8 sites in fig.|^. We 
can see that Uc is finite between 8 = and — 7r/2, taking 
its maximum at Ae ~ 0.5 (9 ~ — 7r/4), or, equivalently, 
Ae 0.5 with tx,t 2 < 0. 

Finally let us make a brief comment. If we change 
t-i from positive to negative, the flat band becomes the 
lowest band, where the ferromagnetic ground state is re- 
alized for arbitrary U > when the flat band is half-filled 
as in Tasaki's model or Mielke's model. 

To summarize, we have proposed that the flat-band fer- 
romagnetism can occur in a non-frustrated structure that 
can model atomic-scale quantum wires. Further investi- 
gations with the density matrix renormalization group 
method will be published elsewhere. Also, the fabrica- 
tion of the As/Si(100) wire is experimentally under way. 
R.A. is grateful to K. Kusakabe for illuminating discus- 
sions. Numerical computations were done on FACOM 
VPP 500/40 at the Supercomputer Center, Institute of 
Solid State Physics, University of Tokyo. 
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FIG. 3. The difference 

in energy, A = £ m in(Smax) — £ : min(0) for 6 electrons in 8 
sites (squares), and for 12 electrons in 16 sites (circles). In 
the inset, we plot the difference in energy, -B m i n (l) — £? min (0) 
(squares) and _E m j n (2) — -E m i n (0) (circles) for 12 electrons in 
16 sites. 



FIG. 1. A tight bindig model proposed by Yajima et al. 
We have indicated the eigenstate on the flat band. Open cir- 
cles represent As atoms, while filled circles Si. Shaded area 
depicts a 'Wannier' state. 



FIG. 4. A similar plot as in Fig. 3 for 

A = -E m in(S' m ax) — -Bmin(O) when an extra transfer 
5t(= O.Fcircles or 0.2:squares) is introduced (as in the inset) 
for 6 electrons in 8 sites. 



FIG. 2. The dispersion relation for the tight-binding FIG. 5. Uc as a function of Ae for 6 electrons in 8 sites, 
model on the connected square lattice for a general value of 
9 (a) and for 9 = -tt/6 (Ae = 0) (b). 
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